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Abstract. New constraints on stellar models are provided by large surveys of massive stars, 
(^~>) ■ interferometric observations and asteroseismology. After a review of the main results so far 

■ obtained, we present new results from rotating models and discuss comparisons with observed 
, features. We conclude that rotation is a key feature of massive star physics. 
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. 1. Large surveys of massive stars 

f~i , Most of the developments in stellar physics arise from the necessity to better reproduce 
O^' observed features. For instance the observation of strong nitrogen enrichments at the 
Q . surface of main sequence OB stars indicates that some extra mixing process is at work. 
^ Since massive stars are fast rotators and since rotation triggers many instabilities able to 
, drive the transport of chemical species, a lot of effort was put in order to properly modelize 
the effects of rotation. Results of these computations at their turn triggered large surveys 
of massive stars made possible thanks to the advent of powerful multispectrographs. 
Among the most recent surveys let us mention the following ones: 

• Keller (2004) presents measurements of the projected rotational velocities of a sample 
of 100 early B-type main-sequence stars in the Large Magellanic Cloud (LMC). The 
sample is drawn from two sources: from the vicinity of the main-sequence turnoff of 

P\J young clusters (ages 1-3 • 10^ y) and from the general field. 

, • Strom et al. (2005) measured the projected rotational velocities for 216 B0-B9 stars 
in the rich, dense h and x Persei double cluster and compared with the distribution of 
OO rotational velocities for a sample of field stars having comparable ages ( 12-15 Myr) and 
O masses (M 4-15 M©). 

J> , • Huang & Gies (2006a) present projected rotational velocities for a total of 496 OB 
' stars belonging to young galactic clusters. Surface helium abundances are given in Huang 
rS & Gies (2006b). 

d • Martayan et al. (2006; 2007b) have measurements of v sin i of B and Be stars in the 

SMC (202 B, 131 Be) and in the LMC (121 B, 47 Be). 

• Wolff et al. (2007) have measured projected rotational velocities for about 270 B-type 
stars belonging to galactic clusters or associations. 

• The VLT-Flames survey (still in progress) has for aim to analyze about 750 OB 
stars observed in 7 fields centered on young clusters, 3 in the Milky Way (NGC 3293, 
4755, 6611), 2 in the LMC (NGC 2004, Nil) and 2 in the SMC ( NGC 330, NGC 346). 
The spectral classification and the radial velocities are discussed in Evans et al. (2005; 
2006) together with some considerations on the populations of binaries and of Be stars. 
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Figure 1. Left panel: Points witii error bars are observations of SMC stars in tlie direction 
of NGC 346 by Hunter et al. (2007) and of NGC 330 by Trundle et al. (2007). Only stars 
with estimated masses between 9 and 13 Mq have been considered. Continuous and dotted 
lines are the tracks computed by Ekstrom et al. (2008). Long-dashed lines are tracks computed 
by Ekstrom et al. (in preparation, see text). Right panel: Observed points are for stars with 
estimated masses between 10 and I6M0 taken from the same references as in the left panel. 
Long-dashed lines are rotating tracks computed for the present work with initial abundances 
as given by Hunter et al. (2007) (see text). The continuous lines are models with rotation and 
magnetic fields (an H follows the initial equatorial velocity in that cases). 

Rotation velocities are discussed in Mokiem et al. (2006; 2007), Dufton et al. (2006) 
and Hunter et al. (2008a) (609 measurements published so far). First results concerning 
the surface chemical compositions are presented in Hunter et al. (2007; 2008b), Trundle 
et al. (2007) (110 measurements of CNO and Mg Si published so far). Helium surface 
enrichments at the surface of O and early B-type stars are discussed in Mokiem et al. 
(2006; 2007) (57 measurements published so far). 

Concerning the distribution of rotational velocities, the following results have been 
obtained: 

• O-type stars in the SMC: Mokiem et al. (2006) deduce the underlying v dis- 
tribution of the uncvolvcd SMC O-type stars. They obtain a mean velocity v of about 
150-180 km s^^ and an effective half width of roughly 100-150 km s^^. 

• OB-type stars in the MCs: Hunter et al. (2008a) present the atmospheric param- 
eters and the projected rotational velocities for approximately 400 O- and early B-type 
stars in the Magellanic Clouds. The observed wsini distributions can be modeled by 
Gaussians with a peak at respectively 100 and 175 km s~^ for the LMC and SMC and 
with a 1/e half width of 150 km in both cases. 

• B- and Be-type stars in the MCs: Martayan et al. (2006; 2007a) obtain mean 
wsini of 161±20 km s^^ and 155±20 km s~^ for SMC B-type stars of respectively 2-5 
(111 stars) and 5-10 M0 (81 stars). Analogous stars in the LMC have mean projected 
rotational velocities of 144±13 and 119±11 km s^^. Let us note that the average velocities 
of Be stars are much greater (see the details in the above references) . As an illustrative 
example, the mean wsini for SMC Be stars in the mass ranges 2-5, 5-10, 10-12 and 
I2-I8M0 are respectively 277±34 (14 stars), 297±25 (81 stars), 335±20 (13 stars) and 
336±40 (14 stars). 

• B-type stars in the Galaxy: Huang & Gies (2006a) present projected rotational 
velocities for 496 OB stars belonging to 19 young galactic clusters with estimated ages 
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between 6 and 73 Myr. Mean vsmi values of 139, 154 and 151 km s^^ have been obtained 
for groups of 09.5-B1.5, B1.5-B5.0 and B5.0-B9.0 type stars. These authors derived the 
underlying probability distribution for the equatorial velocities v and obtained a peak at 
200 km s~^. Dufton et al. (2006) obtain a peak of v at 250 km s~^ with a full- width- half- 
maximum of approximately 180 km s~^ for the unevolved targets in the galactic clusters 
NGC 3293 and 4755. 

Some authors have discussed the variation of the velocity as a function of the age of the 
stars. The main results are the following: 

• O-type stars in the SMC: According to Mokiem et al. (2006) who analyses O- 
type stars in the SMC, the observed distribution of vsini for evolved stars (luminosity 
classes I-II) contains relatively fewer fast rotators and slow rotators compared to the 
distribution for unevolved stars (luminosity classes IV- V). A similar trend is obtained 
for galactic stars. These authors suggest that when the star evolves, it undergoes a spin 
down due to an increased radius, a loss of angular momentum through stellar winds. 
This may explain the smaller proportion of fast rotators. The smaller proportion of slow 
rotators might be due to excess of turbulent broadening among evolved stars. 

• OB-type stars in the MCs: Hunter et al. (2008a) from a sample of 400 O- and 
early B-type stars in the Magellanic Clouds also find that supergiants are the slowest 
rotators in the sample, typically having rotational velocities less than 80 km s~^. 

• OB-type stars in the Galaxy: Huang & Gics (2006b) show that all OB stars of 
their sample experience a spin-down during the MS phase. A few relatively fast rotators 
are found near the TAMS. According to these authors, these stars may be spun up by a 
short contraction phase or by mass transfer in a close binary. 

• Wolff et al. (2007) find that independent of environment, the rotation rates for stars 
in the mass range 6-12 Mq do not change by more than 0.1 dex over ages between about 
1 and 15 Myr. 

Rotation may also depend on the mass and on the metallicity as is deduced from the 

results below: 

• Higher masses, lower velocities: According to Hunter et al. (2008a) there is 
some evidences that the most massive objects rotate slower than their less massive coun- 
terparts. DTifton et al. (2006) find that the mean rotational velocity of stars which have 
strong winds is lower than that of the lower mass stars. 

• Lower Z, higher velocities: Mokiem et al. (2006) find that among O-type stars, 
the distribution for unevolved SMC objects shows a relative excess of stars rotating 
with projected velocities between 120 and 190 km s^^ compared to analogous velocity 
distributions in the Galaxy. This can be interpreted as a decrease of angular momentum 
loss by stellar winds in lower metallicity environments. Hunter et al. (2008a) also obtain 
that SMC metallicity stars rotate on average faster than galactic ones (mainly field 
objects). No difference is found between galactic and LMC stars. Martayan et al. (2007a) 
find that, for B and Be stars, the lower the metallicity, the higher the rotational velocities. 

Many authors find that rotational v(^locities of OB stars in clusters arc greater than 
those of stars belonging to less dense systems like stellar associations or the field: for 
instance 

• Keller (2004) obtains that the mean v sini of early B-type stars in clusters with ages 
between 1-3-10'' years of the LMC is 146 km s~^, while it is 112 km s~^ for analogous 
field stars. A same trend has been found for galactic stars although with lower values for 
both the clusters (116 km s~^) and the field (85 km s~^). 
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Figure 2. Left panel: Observed points are from Hunter ct al. (2007). Tiie long-daslied tracks 
arc from tlie rotating models of Ekstom et al. (in preparation, see text) with Uini/ucrit = 0.4, 
the dot-short dashed tracks are rotating models with Vini = 300 km s~^ from Maeder & Meynet 
(2001). The models of Ekstom et al. have initial values of fiini/ficrit = 0.57, while the models 
of Maeder & Meynet (2001) have nini/Qcrit = 0.64, 0.61 and 0.58 for respectively the 12, 15 
and 2OM0. Right panel: Points with error bars are observations of SMC stars in the direction 
of NGC 346 by Hunter et al. (2007) and of NGC 330 by Trundle et al. (2007). Only stars with 
estimated masses between 15 and 21 M0 have been considered. The star NGC330-003 has not 
been plotted being well outside the range of the other stars, it would lie at the position given by 
log Tcft =4.235 and Alog N/H=1.19. Contirmous and dotted lines are the tracks computed by 
Ekstrom et al. (2008). Long-dashed lines are tracks computed by Ekstom et al. (in preparation, 
see text). 

• Strom et al. (2005) find that B-type stars members of h and x Per (age between 12 
and 15 Myr) have mean usini higher than analogous field stars. The difiPerence between 
these two means depends on the evolutionary stage. For less evolved stars (4-5 Mq), the 
mean projected velocity is 183 km s^^ for cluster stars and 92 km s^^ for field stars, for 
somewhat more evolved stars (5-9 M0), the cluster and field mean vsini are 145 and 93 
km s~^, while for stars approaching the end of the Main-Sequence phase (9-15 Mq), one 
has respectively 104 and 83 km s^^. 

• Dufton et al. (2006) find that the projected velocities in the galactic clusters NGC 
3293 and 4755 are systematically larger than those for the field. Huang & Gies (2006a) 
find from their study of galactic stars that there are more fast rotators among the B 
cluster stars than in the case of the field stars. The mean projected rotational velocities 
are 148±4 km s~^ and 113±3 km s~^ for the cluster and field stars respectively. 

• Wolff et al. (2007) obtain that stars formed in high-density regions lack the cohort 
of slow rotators that dominate the low-density regions and young field stars. 

The Be stars are stars surrounded by an expanding equatorial disks problably produced 
by the concomitant effects of both fast rotation and pulsation. These objects are won- 
derful laboratories to study the effects of extreme rotation. Many new results concerning 
them have been obtained: 

• Objects with Be phenomena are the fastest rotators in the sample studied by Hunter 
et al. (2008a) (400 OB stars in the MCs). This trend is confirmed by Martayan et al. 
(2006; 2007a) who obtain that Be stars rotate faster than B stars whatever the metallicty 
is (see above). 

• These last authors obtain that Be stars with masses below about I2M0 are mainly 
observed in the second part of the MS whatever the metallicity. The more massive stars 
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are mainly in the first part of the MS in the MW, while in the Magellanic Clouds, they 
are all in the second part of the MS. 

• Martayan et al. (2007b) find 13 Be stars among the sample of Be SMC stars with 
short-term periodicity and 9 of them are multi-periodic pulsators. The detected periods 
fall in the range of slowly pulsating B-type stars modes (from 0.4 to 1.60 days). 

• Maeder et al. (1999) and Wisniewski & Bjorkman (2006) find that the fraction of 
Be stars with respect to the total number of B and Be stars in clusters with ages (in 
years) between 7.0 and 7.4 (in logarithm) increases when the metallicity decreases. This 
fraction passes from about 10% at solar metallicity to about 35% at the SMC metallicity. 
These significant proportions imply that Be stars may drastically affect the mean v sin i 
obtained for a given population of B-type stars. 

Surface enrichments in helium and nitrogen have been observed with the following 
main trends (sec also Maeder in these proceedings): 

• He in O-type stars in the SMC: In the SMC, for 31 0-type stars, Mokiem et al. 

(2006) find values of y = ?^Ho/ (n-H + «He) between 0.09 and 0.24, where nne is the density 
number of helium and nu of hydrogen . Note that nuc/in-u + tt-hc) = Y/(Y + iX), where 
Y and X are respectively the mass fraction of helium and of hydrogen. Setting Y ^ 1—X, 
one obtains values for Y = Ay/{l-\-2,y) between 0.28 and 0.56 (here 0.28 would correspond 
to y=0.09 i.e. to the initial helium mass fraction). These authors conclude that while 
rotation can qualitatively account for such enrichments, the observed enrichments are in 
many cases much stronger than those predicted by the models. 

• He in O-type stars in the LMC: In the LMC, for 28 O-type stars, y values 
between about 0.09 and 0.28 are obtained by Mokiem et al. (2007), i.e. helium mass 
fractions between 0.28 and 0.61. 

• He in OB-type stars in the MW: Huang & Gies (2006b) determine He abun- 
dances for OB galactic stars. In their high mass range (8.5Mo < M < I6M0), the He 
enrichment progresses through the main sequence and is greater among the faster ro- 
tatorsf^. On average He abundance increases of 23% ±13% between ZAMS and TAMS. 
These authors also obtain that He enrichments are higher for higher wsini values. Lyu- 
bimkov et al. (2004) find a ZAMS to TAMS increase in He abundance of 26% for stars 
in the mass range 4-11M0 and 67% for more massive stars in the range 12-19M0. 

• N in O-type stars in the SMC: Heap et al. (2006) study a sample of 18 O-type 
stars in the direction of the SMC cluster NGC 346. The surface of about 80% of the stars 
is moderately to strongly enriched in nitrogen, while showing the original helium, carbon 
and oxygen abundances. 

• N in OB-type stars in the MCs: At the present time, the published values of 
nitrogen surface abundances from the VLT-Flames survey concern stars with relatively 
low V sin i (see Trundle et al. 2007 and Hunter et al. 200713) ■ Some of these observations 
for SMC and LMC stars are shown in Figs. [T] and [S] A large spread of abundances is 
found spanning a range between and 1.19 dex of N/H enhancements for masses between 
9 and 21Mq. 

• N in B-type stars in the MW: The galactic stars do not seem to present large 
spread of nitrogen abundances like that seen for Magellanic Cloud stars. Trundle et al. 

(2007) note that if the galactic stars underwent the same degree of enrichment as the 
LMC and SMC stars (in absolute value), this would amount to only a factor of two 

f These authors also found many helium peculiar stars (He- weak and He-strong). These stars 
were not used to study the process of He-enrichment. 

X In Hunter et al. (2008b) nitrogen surface abundances for stars with high v sin i are discussed 
but no detailed tables with the individual measurements are provided. 
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Figure 3. Left panel: Variation of the number ratios of WN to WC stars as a function of 
metallicity. The grey area encompasses the observed ratios. Individual measures are indicated 

by black circles labeled with the name of the Galaxy (see references in Mcynct & Maedcr 2005). 
Solar (0/H) value is taken from Aspluud ct al. (2005), the (0/H) values for the SMC and 
LMC are taken from Hunter et al. (2007). The dotted lines show the predictions of the rotating 
and non-rotating stellar models of Mcynct & Maedcr (2005). The black pentagon shows the 
ratio predicted by Z=0.040 models computed with the metallicity dependence of the mass loss 
rates during the WR phase. The open triangle shows the WC/WN ratio obtained from the new 
rotating models (see text). Right panel: Continuous lines show the evolution of the nitrogen 
surface enrichment as a function of the helium surface enrichment. Model 1 is the fast rotating 
60M© model computed by Meynet & Maeder (2007) (Z = 0.002, Vi^i =523 km s"^ This model 
follows a nearly "homogeneous" evolution. Model 2 is a 20Mo model for Z = 0.002 with an 
initial velocity on the ZAMS of 304 km s"\ Models 3 and 4 are 13M0 for Z = 0.005 with 
an initial velocity on the ZAMS of 382 km s^^ without and with a magnetic field. The dotted 
line shows the typical evolution at the centre for the initial composition used to compute the 
Z = 0.002 models, the dashed line shows the evolution at the centre when the initial composition 
of the LMC as given by Hunter et al. (2007) are adopted. For model 1 the whole track shown 
occurs during the MS phase, for model 2, the star indicates the end of the MS phase, for tracks 
3 and 4 the track ends at the end of the MS phase. The grey area shows the He enrichment 
obtained by Huang & Gies (2006b) at the end of the MS. The three dots correspond to three 
stars observed in the direction of the LMC cluster Nil by Mokiem et al. (2007) and Hunter et 
al. (2007) for which He enrichments have been found Prom top to bottom, one has the stars 
Nll-008 (BO.TIa, ~ 3IM0), Nll-072 (B0.2ni, ~ I4M0) and Nll-042 (BOlll, ~ 25M0). 



or 0.3 dex enhancement in the Galaxy. Such enhancements are similar to the degree of 
uncertainties in their measurements. 

The problem of the mass discrepancy (see Herrero et al 1992), i.e. of the difference 
obtained between spectroscopic and evolutionary masses, has been significantly reduced 
thanks to improvements brought to stellar atmosphere models. However still some dis- 
crepancies are reported in general linked with strong helium surface enrichments: 

• In SMC: Mokiem et al. (2006) find a mild mass discrepancy for stars with spec- 
troscopic masses inferior to about 2OM0, which correlates with the surface helium abun- 
dance. These authors find that the discrepancies are consistent with the predictions of 
chemically homogeneous evolution. Most of the stars observed by Heap et al. (2006) 
exhibit the mass discrepancy problem although no surface He enrichment. 

• In LMC: Mokiem et al. (2007) from the analysis of 0-type stars in the LMC find 
that bright giants and supergiants do not show any mass discrepancy, regardless of the 
surface helium abundance. In contrast they find that the spectroscopically determined 
masses of the dwarfs and giants are systematically smaller than those derived from non- 
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rotating evolutionary tracks. All dwarfs and giants having y > 0.11 {Y > 0.33) show this 
mass discrepancy. 

Close binary stars have a very different evolution from single stars. Some clues pointing 

to these differences are indicated below: 

• Rotational velocities: Hunter et al. (2008a) find that binaries tend to rotate slower 
than single objects. This results however may be somehow biased by the fact that binarity 
is easier to detect at low projected rotational velocities. Huang & Gics (2006a) on the 
contrary find that their binary candidates {i.e. those stars with a difference of radial 
velocity larger than 50 km s~^), have a mean vsini of 173±15 km s~^, higher than 
144±5 km s"'^ the mean for the remaining constant radial-velocity stars. Huang & Gics 
(2006b) find that close binaries generally experience a significant spin-down around the 
stage where polar gravity is equal to 3.9 in logarithm. According to these authors, that 
is probably the result of tidal interaction and orbital synchronization. 

• Surface abundances: According to Hunter et al. (2007) main-sequence binary 
objects have close to baseline nitrogen surface abundances. These systems thus do not 
present apparent signs of extra-mixing. In contrast several evolved binary objects have 
high nitrogen enhancements. These abundances are similar to those observed in appar- 
ently single stars. Thus it appears difficult to discriminate among the possible causes of 
these enrichments, i.e. between extra-mixing operating in single stars and mass-transfer 
events in close binary systems. The same result has been obtained by Trundle et al. 
(2007). 

The above listed results are of course not exhaustive and many more might have 
been cited as the variation of the effective temperature scale for OB stars at different 
metallicities or the variation of the mass loss with the metallicity. Concerning this last 
point let us just mention that Mokiem et al. (2006) find that for stars with logL/L© 
superior to about 5.4, the wind strengths are in excellent agreement with the theoretical 
predictions of Vink et al. (2001). 

Obviously the evolution of rotation and of the surface enrichments depend on the mass, 
the age, the metallicity, the environment (field/clusters), binarity and probably on other 
factors as magnetic fields. Thus the task to disentangle observationally all these different 
effects is very challenging and requires in addition to large surveys detailed observations 
of a few systems for which many precise measurements can be performed. 

2. Interferometry and asteroseismology 

In addition to the above large surveys, there are at least two observational techniques 
which now begin to be applied to massive stars. The first one is interferometry. This 
technique allows to measure the shape of stars, the variation with the colatitude of the 
effective temperature, the inclination angle, the shape of the stellar wind as well as some 
characteristics of stellar disks. Among recent very interesting results let us mention: 

• Meilland et al. (2007b) present the global geometry of the disk around the Be star a 
Arae. The global geometry of the disk is compatible with a thin keplerian disk and polar 
enhanced winds (see also Kervella & Domiciano de Souza 2006). They also obtain that 
a Arae is rotating very close to its critical rotation (for results on the Be star k CMa see 
Meilland et al. 2007a). 

• Domiciano de Souza et al. (2007) have performed the first high spatial and spectral 
resolution observations of the circumstellar envelope of a B[e] supergiant (CPD-572874). 

• Millour et al. (2007), using AMBER/ VLTI observations of the Wolf-Rayet and O 
(WR+0) star binary system 72 Velorum deduce that the binary system lies at a distance 
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Figure 4. Left panel: Variation of tlie mass flux at tlie surface of an initial 35M0 at a stage during 
the core H-burning pliase wlien tlie mass fraction of hydrogen at the centre is Xc — 0.42. The axis 
are in units of solar radius. The velocity of the star on the ZAMS is 550 km s~^ corresponding to 
f2/r2crit ~ 0.84. The star follows a homogeneous evolution. At the stage represented Jl/ficrit ~ 1. 
Right panel: Same as the left panel, but for a later stage with Xc = 0.02 and Sl/ficrit ~ 1. 
Courtesy of C. Georgy. 

of 368+38-13 pc, in agreement with recent spectrophotometric estimates, but significantly 
larger than the Hipparcos value of 258+41-31 pc. 

• Weigelt et al. (2007) have made the first NIR spectro-interferometry of the LBV rj 
Carinae. Their observations support theoretical models of anisotropic winds from fast- 
rotating, luminous hot stars with enhanced high- velocity mass loss near the polar regions. 

The second technique, the asteroseismology, provides new insights on massive star 
interiors (see the paper by C. Aerts in these proceedings). At the moment, data for five 
B-type stars with masses between 8 and 14 M0 have been obtained. The core overshoot 
parameter expressed in units of pressure scale height has been found to be of the order 
of 0.20 (two stars are compatible with that value, two with 0.10 and one with 0.44). For 
three stars the values of the ratio of core to envelope angular velocity have been obtained. 
The values are 5, 3.6 and 1 (solid body rotation). 

In the next sections we discuss a few results recently obtained from massive star ro- 
tating models. 

3. The WC/WN number ratios 

The variation with metallicity of the number of WC to WN stars has been often 
discussed in this meeting. It is a well known fact that this ratio increases with the 
metallicity (see Fig. |3] left paneQ). Many attempts have been performed to reproduce the 
observed trend: for instance the enhanced mass loss rate models of Meynet & Maeder 
(1994) provided a good agreement for solar and higher than solar metallicity but produced 
too few WN stars in metal-poor regions. The inclusion of rotation together with reduced 
mass loss rates accounting for the effects of clumping improved the situation in the metal 
poor region, but produced too many WN stars at solar and higher metallicities (Meynet 
& Maeder 2005). Eldridge & Vink (2006) show that models that include the mass- loss 
metallicity scaling during the WR phase closely reproduce the observed decrease of the 
relative population of WC over WN stars at low metallicities. However such models 
severely underestimate the fraction of WR to 0-type stars. In that case, to improve the 
situation, a high proportion of Wolf-Rayet stars originating from mass transfer in close 
binaries have to be assumed at all metallicities. For instance at solar metallicity about 

f we consider here regions having reached a stationary situation, i.e. regions where the star 
formation rate can be considered to have remained constant for the last twenty million years. 
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Figure 5. Left panel: Evolutionary traclis in tiie HR diagram for a rotating (continuous line) 
and a non-rotating 60Mq stellar model (dashed line) at Z = 0.002. The rotating model has an 
initial velocity of 530 km s~^ or n/ficrit = 0.75. Points indicate the position of the model when 
various values of the hydrogen mass fraction at the centre are reached. Rigtit panel: Evolution 
as a function of the age of the ionising luminosity. Continuous lines are for the rotating model, 
dashed lines are for the non-rotating ones. Courtesy of C. Beffa. 

75% of the WR stars should be produced in close binary systems (see Fig. 5 in Eldridge 
et al. 2008). 

Recently we reexamined this question starting from our rotating stellar models (Meynet 
et al. in preparation). First let us recall that in Meynet & Maeder (2005, 2006), the most 
massive rotating stars enter into the WR regime already during the MS phase. This 
feature has good and bad effects. On one hand, it allows these models to well reproduce 
the variation of the number fraction of WR to 0-type stars since it significantly increases 
the WR lifetimes. On the other hand, it produces very long WN phases since the star 
enters into the WR phase having still a huge H-rich envelope. As a consequence, too low 
values for the WC/WN ratio are obtained at solar and higher metallicities. 

In the above computations, we made the hypothesis that when a star enters into the 
WR stage during the MS phase, it avoids the Luminous Blue Variable phase. This is 
probably not correct. A more realistic solution is to consider that a star which becomes 
a WR star during the MS phase, enters a LBV phase after the core H-burning phase, 
before evolving back into the WR regime. When this more realistic solution is applied, 
reasonable values for both the WR/0 and the WC/WN ratios are obtained. Indeed 
the ratios of WR/0 and of WC/WN given by these models at the solar metallicity 
are 0.06 and 0.9 which compare reasonably well with the observed values of 0.1 and 
0.9 respectively. Both ratios are not reproduced by the non-rotating models to which a 
similar solution is applied. At the moment only the case at solar metallicity has been 
computed, but we are confident that such a scenario will also provide reasonable answers 
at other than solar metallicities. This discussion illustrates the possible key role that the 
LBV phase may play in shaping the WC/WN ratio. 



4. The changes of the surface abundances 

Models without rotation predict no surface enrichments before the red supergiant 
stage for stars less massive than about 4OM0. This is clearly in contradiction with the 
observations. In contrast, rotationally induced mixing produce changes of the surface 
abundances already during the MS phase. This is shown in Figs. [T] and [21 Model results 
depend on the physics included and on the choice of some parameters. In Figs. [1] left and 
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[fright, two series of models are plotted: the models of Ekstrom et al. (2008) have been 
computed with the same physics as in Meynet & Maeder (2005), i.e. the expression for 
the horizontal turbulence was taken as proposed in Zahn (1992), the value of a in the 
expression of the shear diffusion coefficient was taken equal to 2 (see Eq. 3 in Maeder 
& Meynet 2001), and an overshooting parameter of O.lHp was used (these models will 
be called models A hereafter), the models B by Ekstrom et al. (in preparation) differ in 
many aspects. We just mention the two most important ones: they use the expression 
for the horizontal turbulence by Maeder (2003) and the value of a in the expression of 
the shear diffusion coefficient was taken equal to 4. We can clearly see that, given an 
initial value of f2ini/ricrit, models B are more efficiently mixed than models A. One sees 
also that models A, even spanning the whole range of possible fiini/r^crit values cannot 
reproduce the highest enrichments. Models B in contrast do appear in better position to 
reproduce the observed range of values and are thus to be preferred. 

In Fig. [T] right, rotating models computed with the same physics as models B have 
been computed for the LMC metallicity (long-dashed lines). Tracks computed with the 
effects of a magnetic field as in Meynet & Maeder (2007) are also shown for a few 
velocities (continuous lines). We see that models with magnetic fields are more efficiently 
mixed, they also show a more progressive surface enrichment. Both series of models can 
reproduce the observed range of values. One notes however that models with magnetic 
fields appear in a better position for explaining the most extreme cases. In Fig.[2]left, the 
models B (long-dashed lines) are plotted together with the models of Maeder & Meynet 
(2001). Again models B are more efficiently mixed and appear in better position for 
reproducing the observed enrichments. 

On the whole, we see that the recent observations obtained in the large surveys de- 
scribed above support quite efficient mixing processes. The models of Ekstrom et al. (in 
preparation) can account for a great part of the observed range. More detailed compar- 
isons will be presented elsewhere. 

In Fig. [3] right, the continuous lines show the evolutions of the surface N-enrichment 
as a function of the He-enrichment at the surface of different models. The dotted line 
shows the evolution in the convective core of a 2OM0 stellar model with Z = 0.002. The 
evolution in the convective core represents an upper envelope of what we can expect from 
a theoretical point of view. Note that this curve does not much depend on the initial 
mass, but depends on the initial CNO content as can be seen by comparing the dotted 
and the dashed line. The dashed line is obtained from models with initial abundances 
for the SMC as given by Hunter et al. (2007), while the dotted one is obtained starting 
with the relative abundances as given by Asplund et al. (2005). 

Let us emphasize that progression along the "surface" (continuous lines) and "core" 
(dotted) lines goes at a different pace. This can be realized by noting that on the "surface" 
line the end of the MS phase in general occurs very early (typically the end of the MS is 
indicated by a star on the track 2, corresponding to the same model as the one used to 
draw the dotted line), while, on the "core" line, the end of the MS phase would occurred 
at Kj — Yc = 0.75 and Alog N/H tending toward infinity since the hydrogen abundance is 
zero. The N-enrichment occurs very rapidly at the surface well before any surface helium 
enrichment. This is of course due to the fact that very rapidly (see the dotted line) 
the nitrogen abundance increases at the centre creating thus a strong chemical gradient 
between the core and the envelope. Since diffusive velocity is greater when the gradient 
of abundance is greater, the presence of such a strong gradient favors a rapid mixing. 
The gradient in helium is much more shallow and makes the diffusion of this element to 
occur on much longer timescales. 

With these remarks in mind, let us now discuss each model: model 1 shows the behavior 
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of a fast rotating 60 Mq model at Z = 0.002 following a nearly homogeneous evolution. 
The evolution of its surface abundances approaches that of the central ones. The whole 
portion of the track shown in Fig. [3] right occurs during the MS phase. Such a model 
would easily account for the observed He enrichments during the MS phase. The model 
2 (20 Mq with Uini = 304 km s^^ and for Z = 0.002) shows at the end of the MS 
phase an increase in helium of 0.01 (in mass fraction) and an increase of 0.7 dex in N/H. 
The He-enrichment is lower than the observed one. Models 3 and 4 compare two 13 Mq 
models for the LMC composition. Both models have very high initial rotation (382 km 
s~^) but model 3 is computed without magnetic field while model 4 is computed with a 
magnetic field. We see that despite the high velocity, the model without magnetic field 
does not have any surface helium enrichment. Model with magnetic field reaches at the 
end of the MS phase an enrichment of nearly 0.015 with respect to its initial value. This 
enrichment is however too low to account for the He-enrichments obtained by Huang & 
Gies (2006b) at the end of the MS phase. 

We conclude from these comparisons that the present-day estimates of He-enrichments 
at the surface of some stars require tracks following a nearly homogeneous evolutiorlj. 



5. Wind anisotropies and homogeneous evolution 

Rotation has a deep impact on the way massive stars are losing mass (see the papers 
by A. Maeder and R. Hirschi in these proceedings). Here we focus on the anisotropies of 
the winds induced by fast rotation (Maeder 1999). When the surface velocity is near the 
critical one strong wind anisotropies appear. This is illustrated in Figs. |4] which shows 
the variation of the mass flux at the surface of a fast rotating 35 Mq model following 
a homogeneous evolution. The lengths of the arrows are proportional to the mass flux. 
Figure [4] left shows a stage when the star is an 0-type star, the right panel when it is 
a WR star. The scale for the mass flux was changed between left and right panels since 
during the WR phase the mass fluxes are about 2 orders of magnitudes higher. The grey 
scale varies as a function of the effective temperature (von Zeiple effect). 

Accounting for wind anisotropies allows to keep in the star about 25-30% of the angular 
momentum which would be lost by an isotropic wind. This difference produces important 
differences in the way the angular momentum is distributed in the star at the pre- 
supernova stage. Typically, when wind anisotropies are considered, the specific angular 
momentum in the 3 inner solar masses are more than a factor 6 higher than when the 
isotropic winds are considered (see Meynet & Maeder 2007) . Thus wind anisotropies are 
probably a key feature of homogeneous evolutionary tracks and of GRB progenitors. 

Such homogeneous evolution will also produce higher outputs of ionising photons as 
can be seen in Fig. O The total number of photons released during the whole stellar 
lifetime with an energy sufficient to ionize H is about a factor 2 greater (passing from 
1.46 10^^ to 2.65 10^^) for the blueward track than for the normal one. Most of these 
photons (98%) are emitted during the core H-burning phase. The total number of Hel 
ionising photons passes from 2.4 10^^ to 6.7 10^^. Most of these photons are released 
during the core H-burning phase (92%). The number of Hell ionising photons passes 
from 8.7 10^^ to 184 10^^. In that case 2/3 are emitted during the core He burning 
phase. Thus homogeneous evolution increases the number of Hel and HcH photons by 
about a factor 3 and 20 respectively. Therefore, if homogeneous evolution is a not too 

f Note that many stars show N enrichments but no He enrichment, they would stand on the 
vertical line at the abscissa in Fig. [3] right. These stars can in general be explained by models 
with normal rotation velocities. 
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rare scenario for PopIII and very metal poor massive stars, then the budget of ionising 
photons in the early Universe should take into account such sources. 

6. Conclusions 

Considering a massive star of given initial mass and metallicity, one encounters for 
increasing angular momentum content the following effects of rotation: The first effect of 
rotation which already occurs for modest rotation rate is internal mixing. Of course when 
the velocity increases the mixing also increases. Then increasing the angular momentum 
content, the second series of effects induced by rotation concerns mass loss. Rotation may 
trigger mechanical mass loss, it may also at very low metallicity increase the "metallicity" 
of the outer layers (metals produces by the star itself) and thus increase the opacity and 
trigger radiative stellar winds (see the papers by Hirschi and Ekstrom in the present 
volume). Finally as an extreme case of the effects of rotation, there are the homogeneous 
evolution. Theory indicates that many of these effects vary as a function of Z, tending to 
be more important in metal poor regions. At the moment many indirect observational 
features seem to require high rotation rates at low Z for massive stars (see Meynet et al. 
2008). Direct observations, limited to MC metallicities, seem also to support this view 
by indicating that rotation rates appear to be higher at lower Z. Therefore rotation is 
probably a key effect for understanding the first generations of massive stars. 
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Discussion 

KOENIGSBERGER: Could you say a few words on how the gradient of mean molecular 
weight can affect the efficiency of rotational mixing? 

Meynet: When the molecular weight increases with depth as is the case in stars, mixing 
becomes more difficult, since greater energy is required to lift off heavy material and to 
mix down light one. This can be shown through the expression of the Richardson criterion 
(see Maeder & Meynet 1996, A&A, 313, 140). Meynet and Macder (1997, A&A, 321, 165) 
have shown that the strict application of the Richardson criterion would prevent any 
mixing in regions with a molecular weight gradient. Only when the effects of the strong 
horizontal turbulence arc accounted for in a proper way (Talon & Zahn 1997, A&A, 317, 
749) can mixing be efficient enough. 

LiMONGi: A comment on the anticorrclations obs(TV(xi in globular clusters. In principle 
AGB models can explain the anticorrclations maintaining the C+N+0=const. In the 
most massive AGB stars, HBB occurs and also if there are few dredge-up episodes the 
chemical composition observed in the turn off stars of globular clusters can be reproduced, 
the problem is that to quantitatively account for the observations a very peculiar IMF is 
required. Now a question concerning the WC/WN ratio: how much the WN/WC problem 
depends on the mass loss rate during the red supergiant phase? 

Meynet: The WR/0 and the WC/W'N ratios can be reproduced without modifying 
the mass loss rates in the red supergiant phase. Thus, the mass loss rate during the 
RSG does not appear as a critical point. However higher mass losses during the red 
siipergiant phase might allow lower initial masses to enter the WN regime and thus may 
have an impact on this ratio. At the moment however, except may be in the clusters 
near the galactic centre, there is no single-aged clusters showing both a significant red 
supergiant and W^olf-Rayet population. This indicates that the two populations originate 
from different initial mass ranges. 

Hunter: While fast rotation does change the surface gravities, it does not change the age 
of all the stars and hence the bulk of the "older" fast rotating unenriched stars remain. 
Secondly, given the rotational velocity distribution, the low vsmi stars are not due to 
sini effects, but are indeed slowly rotating stars. 

Aerts: 1) The study of Morel et al. (2006) contains stars with N-enrichment and a 
seismic estimate of the equatorial rotational velocity, so here one is not bothered with 
sinz uncertainty. 2) The FLAMES survey contains two clusters with a lot of /3Cephei 
pulsators, so a line-profile analysis of these stars can also lead to sini. this could help to 
interpret the "Hunter-diagram" . 



